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ABSTRACT 

We measure the topology of the main galaxy distribution using the Seventh Data Release of the Sloan 
Digital Sky Survey, examining the dependence of galaxy clustering topology on galaxy properties. The 
observational results are used to test galaxy formation models. A volume-limited sample defined by 
Air < —20.19 enables us to measure the genus curve with amplitude of G = 378 at 6ft-~^Mpc smoothing 
scale, with 4.8% uncertainty including all systcmatics and cosmic variance. The clustering topology 
over the smoothing length interval from 6 to 10/i~"*Mpc reveals a mild scale-dependence for the shift 
{Av) and void abundance (Ay) parameters of the genus curve. We find substantial bias in the topology 
of galaxy clustering with respect to the predicted topology of the matter distribution, which varies 
with luminosity, morphology, color, and the smoothing scale of the density field. The distribution 
of relatively brighter galaxies shows a greater prevalence of isolated clusters and more percolated 
voids. Even though early (late)-type galaxies show topology similar to that of red (blue) galaxies, the 
morphology dependence of topology is not identical to the color dependence. In particular, the void 
abundance parameter Ay depends on morphology more strongly than on color. We test five galaxy 
assignment schemes applied to cosmological N-body simulations of a ACDM universe to generate 
mock galaxies: the Halo-Galaxy one-to-one Correspondence model, the Halo Occupation Distribution 
model, and three implementations of Semi- Analytic Models (SAMs). None of the models reproduces 
all aspects of the observed clustering topology; the deviations vary from one model to another but 
include statistically significant discrepancies in the abundance of isolated voids or isolated clusters and 
the amplitude and overall shift of the genus curve. SAM predictions of the topology color-dependence 
are usually correct in sign but incorrect in magnitude. Our topology tests indicate that, in these 
models, voids should be emptier and more connected, and the threshold for galaxy formation should 
be at lower densities. 
Subject headings: large-scale structure of universe - cosmology: observations - methods: numerical 



L INTRODUCTION 

Galaxy clustering has long been used to constrain cos- 
mological models and to understand formation of galax- 
ies. The most extensively studied clustering statistics 
are the autocorrelation function and the power spectrum. 
They are Fourier transforms of each other, and measure 
the clustering strength as a function of scale. The "Stan- 
dard Cold Dark Matter" model (with the density param- 
eter Qm = 1, scale-invariant primordial fluctuations, and 
standard relativistic particle background), popular in the 
1980s, was ruled out by showing that it was inconsistent 
with the observed galaxy correlation function (Maddox 
et al. 1990) and power spectrum (Vogeley et al. 1992; 
Park et al. 1992, 1994). These statistics are also used to 
measure the biasing in the galaxy clustering amplitude 
with respect to matter, and are key constraints on galaxy 
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formation models connecting dark matter halos and lu- 
minous galaxies such as semi-analytic models (SAM; e.g.. 
White & Frenk 1991; Kauffmann et al. 1993; Cole et al. 
1994, 2000; Benson et al. 2002; Bower et al. 2006; Cat- 
taneo et al. 2006; Croton et al. 2006; De Lucia & Blaizot 
2007; Monaco et al. 2007; Somerville et al. 2008) and 
halo occupation distribution models (e.g., Seljak 2000; 
Peacock & Smith 2000; Scoccimarro et al. 2001; Berlind 
& Weinberg 2002; Kang et al. 2002; Berhnd et al. 2003; 
Zheng, Coil, & Zehavi 2007). 

Topology analysis was introduced by Gott et al. (1986, 
1987) to test the Gaussianity of the primordial density 
fluctuations, which is one of the key characteristics of 
simple inflationary models (Bardeen et al. 1986). At 
large scales density fluctuations are still in the linear 
regime and maintain their initial topology, and it is pos- 
sible to check whether or not the primordial fluctuations 
were a Gaussian field. At smaller scales, non-linear grav- 
itational evolution and biased galaxy formation make 
the topology of the observed galaxy distribution devi- 
ate from the Gaussian form even if the initial conditions 
were Gaussian distributed as shown by perturbation the- 
ories and large N-body simulations (Park & Gott 1991; 
Weinberg & Cole 1992; Matsubara 1994; Matsubara & 
Suto 1996); using fractional volume rather than density 
threshold as the independent variable in topology anal- 
ysis mitigates but does not eliminate these non-linear 
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and biasing effects (Weinberg et al. 1987; Melott et al. 
1988). Through studies of many observational samples, 
the topological properties of the large-scale distribution 
of galaxies have been examined (Gott et al. 1989; Park, 
Gott & da Costa 1992; Moore et al. 1992; Vogeley et 
al. 1994; Rhoads, Gott, & Postman 1994; Protogeros 
& Weinberg 1997; Canavezes et al. 1998; Park, Gott 
& Choi 2001; Hoyle, Vogeley & Gott 2002; Hikage et 
al. 2002, 2003; Park et al. 2005; James, Lewis & Col- 
less 2007; Gott, Choi & Park 2009; James et al. 2009). 
On non-linear or quasi-linear scales, topology analysis is 
useful in constraining both cosmological parameters and 
galaxy formation mechanisms (Park, Kim & Gott 2005; 
Gott et al. 2008). In particular, differences in cluster- 
ing topology for different types of galaxies reflect their 
different history of formation and evolution. Therefore, 
looking at the topology of large-scale structure traced by 
different types of galaxies can put strong constraints on 
galaxy formation mechanisms. 

In this paper we use the Seventh Data Release (DR7; 
Abazajian et al. 2009) of the Sloan Digital Sky Sur- 
vey (York et al. 2000) to measure the topology of the 
galaxy distribution and its dependence on galaxy lumi- 
nosity, color, and morphology. DR7 constitutes the final 
release of the SDSS Legacy Survey, and thus of the main 
SDSS galaxy redshift survey. We supplement the SDSS 
data with missing redshifts to increase the completeness 
of the redshift catalog. We then generate a set of volume- 
limited samples of the SDSS galaxies divided according 
to their luminosity, morphology, and color to study the 
relation between the topology and properties of galaxies 
tracing the large-scale structure. 

2. THE KIAS-VAGC CATALOG 

We prepare a catalog containing 593,514 redshifts of 
the SDSS (York et al. 2002; Stoughton et al. 2002) 
main galaxies (Strauss et al. 2002) in the magnitude 
range of 10 < r < 17.6, which is called the Korea Insti- 
tute for Advanced Study (KIAS) Value-Added Galaxy 
Catalog (VAGC) (Choi et al. 2010). Its main source is 
the New York University VAGC Large Scale Structure 
Sample (brvoidO) from Blanton et al. (2005), which 
lists redshifts of 583,946 galaxies with 10 < r < 17.6 
(8,562 with 10 < r < 14.5). We excluded 929 objects 
that are in error. They are outer parts of large galax- 
ies that are deblended by the automated photometric 
pipeline, blank fields, or stars, etc. Out of the remaining 
583,017 galaxies, 5.7% of galaxies in the NYU VAGC LSS 
do not have spectra due to fiber collision and have red- 
shifts borrowed from the nearest neighbor galaxy. We 
added redshifts of 10,497 galaxies with 10 < r < 17.6 
(1,455 with 10 < r < 14.5) to the KIAS VAGC that are 
missing in the NYU VAGC LSS. An SDSS photomet- 
ric catalog and various existing redshift catalogs such as 
the Updated Zwicky Catalog (UZC; Falco et al. 1999), 
the IRAS Point Source Catalog Redshift Survey (PSCz; 
Saunders et al. 2000), the Third Reference Catalogue 
of Bright Galaxies (RC3; de Vaucouleurs et al. 1991), 
and the Two Degree Field Galaxy Redshift Survey (2dF- 
GRS; CoUess et al. 2001) are used in this step. During 
this supplementation process, we corrected positions of 
galaxies having wrong central positions, separated merg- 
ing objects, and removed spurious objects such as parts 
of big late- type galaxies. The angular survey mask given 



by NYU- VAGC is maintained and covers 2.562 sr of sky 
with angular selection function greater than 0. While 
the SDSS main galaxy sample extends to r = 17.77 over 
most of the sky, we retain the shallower limit r = 17.60 
used during the early part of the survey (see Strauss et 
al. 2002) to maintain a homogeneous depth of the survey 
volume and thus a less complicated outer boundary. Ex- 
tending to r = 17.77 would add roughly 114, 303 galaxies 
(from NYU VAGC LSS fullO sample). 

The resulting KIAS VAGC has an angular selec- 
tion function significantly better than the original NYU 
VAGC LSS sample in high surface density regions. Fig- 
ure [T] shows the distribution of the SDSS galaxies in the 
northern galactic hemisphere. To minimize boundary ef- 
fects and to keep the surface-to- volume ratio low in our 
topology analysis, we do not use the galaxies in three 
narrow stripes of the southern Galactic cap and in the 
small region containing the Hubble Deep Field, and we 
trimmed some regions with narrow angular extent. What 
remains are the galaxies inside the green solid boundary 
shown in Figure [T] The remaining survey region with the 
angular selection function greater than covers 2.33 str. 
Black points in Figure [1] are the galaxies in NYU VAGC 
LSS catalog and red points are those newly added. It 
can be seen that the added galaxies are mostly located 
in high surface density regions. In particular, a large 
number of galaxies are added in the Virgo cluster re- 
gion and near the equator. Correspondingly, we recalcu- 
late the angular selection function of the KIAS VAGC 
within each spherical polygon formed by the adaptive 
tiling algorithm (Blanton et al. 2003) used for the SDSS 
spectroscopy. Choi et al. (2010) show the angular selec- 
tion function before and after supplementation. After the 
supplementation, the survey area (in the northern hemi- 
sphere) having angular selection function greater than 
0.97 increases from 39.8% to 54.3% of the area with the 
selection function greater than 0. 

All atlas images of galaxies in our KIAS VAGC are 
downloaded from Princeton SDSS reduction serveiQ and 
analyzed to measure the seeing- and inclination-corrected 
i-band Petrosian radius, i-band inverse concentration in- 
dex, Cin and A{g — i) color gradient (Park & Choi 2005; 
Choi et al. 2007). These parameters are included in 
KIAS VAGC together with other fundamental photomet- 
ric and spectroscopic parameters supplied by the NYU 
VAGC. 

All galaxies in KIAS VAGC are classified into early 
(elliptical and lenticular) and late (spiral and irregular) 
morphological types. We use the automated classifica- 
tion scheme developed by Park & Choi (2005) using u — r 
color, A(g — i) color gradient, and the inverse concen- 
tration index Qn. Reliability and completeness of this 
classification reach about 90%. To improve the auto- 
mated classification results by visual inspection we se- 
lected 82,323 galaxies that are located in the "trouble 
zones" where the reliability of the automated classifi- 
cation is low. These are the galaxies having neighbors 
at very small separations (for these cin is inaccurate) or 
those classified as blue early types or red late types. Thir- 
teen astronomers made the visual inspection to check the 
morphological type of these galaxies. As a result, the 
morphology of 7% of the inspected galaxies was changed 

^ http://photo.astro.princcton.edu 
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Fig. 1. — SDSS galaxies in the northern galactic hemisphere. Black points are the galaxies in the NYU VAGC LSS catalog and red 
points are those newly added in the KIAS-VAGC. Green solid lines delineate the boundaries of the analysis regions, (rj, A) are the SDSS 
survey coordinates. Lines of constant RA and DEC lines arc shown in black as roughly vertical and horizontal curves, respectively. 



and sonic spurious objects are removed. We kept the 
morphology of 5,956 galaxies in the SDSS LSS-DR4plus 
sample that were already inspected by Choi et al. (2007). 
Figure 2 shows the galaxies of the KIAS-VAGC con- 
tained in ten consecutive wedge slices in the main survey 
area. Each slice has the radius of 450 /i~^Mpc and is 
6 degree thick in declination. The range in right ascen- 
sion is varied according to the survey boundaries shown 
in Figure 1. The right ascension and declination ranges 
and the number of galaxies in each slice are listed in 
Table 1. Galaxies are distinguished according to their 
color: blue is given to those with u — r < 2.0, green is for 
2.0 < u — r < 2.5; and red is for u — r > 2.5. Since we 
are showing galaxies in the apparent magnitude-limited 
sample, galaxies are on average faint and blue at small 
distances and bright and red at large distances. A part 
of the Sloan Great Wall (Gott et al. 2005) is seen in the 
bottom slice (slice j, at r « 25Qh~^ Mpc), and the CfA 
Great Wall (Geller & Huchra 1989) appears in slice e at 
r « 100/i~^ Mpc. A part of the Cosmic Runner (Park 
et al. 2005) shows up in the top slice a. There are also 
numerous prominent superclustcrs, notably a group of 
superclustcrs in slice / and a few rod-shaped ones in slice 
i. In addition to these large scale overdense structures, 
many large voids are prevalent in the survey volume. It 
should be also noted that none of these features really 
stands out any more; the Sloan Great Wall is just a high 
end of the distribution of structures that are quite com- 
mon. A weak circular symmetry around the observer is 



seen in slice j, but we do not have such impression in 
other slices. 

We will use several volume-limited samples defined by 
r-band absolute magnitude cuts and the redshift limits 
of 2:min = 0.02 and Zmax, where z^ax is determined by 
the survey flux limit r = 17.6 and the absolute mag- 
nitude cut. We frequently refer to the "mean separa- 
tion" between galaxies in a sample, by which we mean 

J = rig where Ug is the number density of galaxies. 
Figure [3] shows the galaxies in our catalog in the redshift 
versus r-band absolute magnitude (Mr) plane and the 
boundary lines defining one of our volume-limited sam- 
ples "BEST" , which contains the maximum number of 
galaxies (133,947) among volume-limited samples. It is 
amusing to note that the BEST sample is about 1000 
times larger than the galaxy sample (drawn from the 
CfAl redshift survey) first used for topological analysis 
by Gott et al. (1986, 1987). 

The WMAP 3 year cosmological parameters of flm = 
0.27, r^A = 0.73, fif, = 0.044, /i = 0.72, n^ = 0.958 
(Spergel et al. 2007) are used to relate redshift and co- 
moving distance. 

3. THE GENUS STATISTIC 

We measure the topology of the galaxy distribution 
using the general methods set out by Gott et al. (1987, 
1989), with statistical characterizations of genus curve 
shape (Eqs. 0] and [S] below) introduced by Park et 
al. (1992, 2005). The point distribution of galaxies 
is smoothed by a constant-size Gaussian filter and iso- 
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Fig. 2. — Distribution of SDSS galaxies in the legacy survey area divided into ten wedge slices. Each diagram has radius of 450 h^^Mpc 
and thickness of 6 degrees in declination. The boundaries of each wedge diagram are defined in Table ??. Galaxies arc color-coded hy u — r 
color. Blue, green, and red dots represent galaxies with u — r < 2.0, 2.0 < u — r < 2.5, and u — r > 2.5, respectively. The radial coordinate 
is the comoving distance in units of h~^Mpc. Right ascension is given in units of hour. A part of the Sloan Great Wall (Gott et al. 2005) 
is seen in slice j. 



density contour surfaces of the smoothed galaxy density 
distribution are searched to calculate the genus. The 
genus is defined as 

G = Number of holes — Number of isolated regions (1) 

in the iso-density contour surfaces at a given threshold 
level, v, which is related to the volume fraction / on the 



high density side of the density contour surface by 
1 



/ 



/27r 



e-^'/^dx. 



(2) 



The / ~ 50% contour corresponds to the median volume 
fraction contour {v = 0). In our analysis we will always 
use the volume fraction to find the contour surfaces and 
the threshold level v is related with the volume fraction 
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the following three variables: 

/_.^g(t/)t/dt/ 



Fig. 3. — Galaxies in tlic KIAS VAGC shown in A'U-z space. The 
green curve is the flux Hmit of r = 17.6, and the straight lines in 
magenta define the sample BEST, which includes the maximum 
number of galaxies among volume-limited samples. Black points 
are the galaxies in the NYU VAGC LSS catalog and red points are 
those newly added in the KIAS- VAGC. Blue points are those to be 
added when the limiting apparent magnitude extends to r = 17.77, 
but they are not used in this study. 



TABLE 1 
Survey Regions for the Wedge Diagrams in Figure [2] 



Panel Name Right Ascension Declination Num. of galaxies 



a 


8.5° - 15.3° 


54° - 60° 


22,237 


b 


8.5° - 16.0° 


48° - 54° 


23,950 


c 


8.0° - 16.0° 


42° - 48° 


27,232 


d 


8.0° - 17.0° 


36° - 42° 


36,962 


e 


8.0° - 17.0° 


30° - 36° 


41,217 


f 


8.0° - 16.5° 


24° - 30° 


47,140 


9 


8.0° - 16.5° 


18° - 24° 


45,728 


h 


8.0° - 16.2° 


12° - 18° 


47,036 


i 


8.8° - 16.0° 


6° - 12° 


45,554 


J 


8.8° - 15.7° 


0° - 6° 


43,088 



Au 



1-2.2 aci'') dv' 



Ar = 



jl'^9G{v)dv' 



(4) 



(5) 



where goiv) is the genus of the best- fit Gaussian curve 
(Eq. 3). Thus Av measures any shift in the central por- 
tion of the genus curve. The Gaussian curve (Eq. 2) has 
/S.V = 0. A negative value of /S.v is frequently called 
a "meatball shift" caused by a greater prominence of 
isolated connected high-density structures that push the 
genus curve to the left (Weinberg et al. 1987). Ay and 
Ac measure the observed number of voids and (super) 
clusters relative to those expected from the best-fitting 
Gaussian curve. 

Because the genus-related statistics are defined as in 
Equations (4) and (5), they depend on the amplitude of 
the best-fit Gaussian curve, which has a statistical uncer- 
tainty. The fluctuation in the amplitude yields numerical 
uncertainties in these statistics. One might want to fix 
the amplitude using the observed power spectrum and 
the relation between the power spectrum and the ampli- 
tude of a Gaussian genus curve (given just after Eq. 3). 
But the power spectrum measured from an observational 
sample also has statistical uncertainties, which propagate 
to the uncertainty in the amplitude of the genus curve. 
Furthermore, the genus versus power spectrum relation 
in the case of Gaussian fields cannot be directly applied 
to non-Gaussian fields in general. We also want to sepa- 
rate the non-Gaussianity seen in the genus curve due to 
the amplitude drop from that due to the deviation of the 
genus curve from the Gaussian shape, and this is best 
done when the observed genus curve is compared with a 
Gaussian curve with their amplitudes matched together. 
When galaxy formation models are tested against obser- 
vations in section 5, we compare the genus- related statis- 
tics calculated in the exactly same way for both obser- 
vational data and the mock samples of galaxy formation 
models. 



by Equation 2. 

For Gaussian random phase initial conditions the genus 
curve is: 

g{v) = A{1 - v')e''^"l\ (3) 
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where the amplitude A = {Ip) jifl"^ J^-k^ and Ip 
the average value of fc^ in the smooth power spectrum 
(Hamilton, Gott, & Weinberg 1986; Doroshkevich 1970). 
Deviation of an observed genus curve from Equation [3] 
indicates the non-Gaussian signal. 

The shape of the genus curve can be parameterized by 
several variables (Park et al. 1992, 2005; Vogeley et al. 
1994). The first is the amplitude of the genus curve as 
measured by the amplitude of the best fitting Gaussian 
curve of Equation [31 This gives information about the 
power spectrum and phase correlation of the density fluc- 
tuation. Deviations in the shape of the genus curve from 
the theoretical random phase case can be quantified by 



4. RESULTS 

4.1. Genus Measurement 

We first study how the topology of large scale struc- 
ture changes as the smoothing scale changes. For this 
purpose we use the volume-limited sample, BEST. It in- 
cludes galaxies with Mr < —20.19 -I- 51og/i (51og/i will 
be subsequently dropped in the expression of Mr) and 
0.02 < z < 0.116, where the upper redshift limit is de- 
termined by the apparent magnitude limit r = 17.6 and 
Mr = —20.19. The mean separation between galaxies in 

BEST is J = rig ^/^ = 6.1/i^iMpc. FigureHshows the co- 
moving number density of galaxies in BEST as a function 
of radial comoving distance R. It can be seen that the 
comoving density is roughly constant, but there is about 
20% excess in the density at i? = 210 -- 250/i"iMpc 
caused by the Sloan Great Wall. A genus curve obtained 
from BEST is shown in Figure [5] and the data are given 
in Table [TOl in Appendix O 
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TABLE 2 
N-BODY Simulation Parameters 
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Fig. 4. — Comoving number density of galaxies in the BEST 
sample as a function of radial comoving distance R. 



To calculate the genus a smooth density field is 
obtained from galaxy positions in comoving space 
smoothed by an Ra = 6.0/i~^Mpc Gaussian, 



'-mh 



This definition differs from 



(27ri?2.)-3/2exp(- 

that of some previous works (Hamilton et al. 1986; Park 

et al. 1992; Vogeley et al. 1994), which is related with 

current definition by A = \/2-Rg- The solid line with 
error bars shows measurements from the observational 
data, and the long-dashed line is its best-fit Gaussian 
curve. As a comparison, we include the results (curves 
with smaller amplitudes) from the BEST sample of SDSS 
Data Release 3 (DR3) measured by Park et al. (2005). 
As the DR7 sample covers the northern sky without a 
gap, the volume-to-surface ratio increased significantly 
since DR3, and the amplitude of the genus curve in- 
creased remarkably. The amplitude is now G = 373 ± 18 
with only 4.6% uncertainty, including all systematic ef- 
fects and cosmic variance. 

Figure [5] shows that the observed genus curve deviates 
from the Gaussian one in such a way that the genus at 
high thresholds has smaller amplitude. This means there 
are fewer voids and fewer superclusters when compared 
with the Gaussian genus curve. In other words, voids 
and superclusters are more connected and their sizes are 
larger than those expected for Gaussian fields. In addi- 
tion, the positive (G > 0) region of the genus curve is 
slightly but systematically shifted towards lower v val- 
ues, and the percolation of large-scale structure (maxi- 
mum G) occurs slightly below z/ = 0. According to the 
perturbation theory prediction of Matsubara (1994), if 
there were only non-linear gravitational effects involved, 
with no biasing, then one should find Ay + Ac = 2 at 
all scales (Park et al. 2005), and the N-body results pre- 
sented for dark matter in Figure fT8l below are consistent 
with this expectation. Since the observed Ay and Ac 
are both less than 1, it shows that biasing effects are 
involved. 

4.2. Mock Surveys and Systematic Biases 

The genus curve in Figure [5] is affected by peculiar ve- 
locity distortions, boundary effects, and finite sampling 



Parameter 




SI 


S2 


N^ 




2048^ 


2048^ 


Nm 




2048^ 


2048^ 


Lbox 




1024/1-1 Mpc 


1433.6/1-1 Mpc 


Nstert 




1880 


1250 


z^ 




47 


50 


m^ 


8.3 


X W^h-^AlQ 


2.3 X lOiO/i-iM© 


7 




0.05/1-1 ^^p^ 


0.07/1-1 Mpc 


2.5 


X lO"/i-iM0 


6.8 X IOII/1-IM0 


dnL 




4.4?t-i Mpc 


5.9/1-1 Mpc 


h 




0.737 


0.737 


Hm 




0.238 


0.238 


fit 




0.042 


0.042 


Ha 




0.762 


0.762 


Us 




0.958 


0.958 


b 




1.314 


1.314 



''Number of CDM particles 

1" Number of global evolution timcsteps from Zi to present epoch 

"^Initial redshift whore the simulation started 

'I Simulation particle mass 

•^ Force resolution 

Minimum halo mass 
^Mean halo separation 



of the density field. However, using detailed comparisons 
to mock catalogs from cosmological simulations, we will 
show that the deviations from the Gaussian field pre- 
diction, and indeed from the expected non-linear matter 
density field, are genuine, and statistically significant. 

We use mock galaxies to estimate uncertainties in the 
measured genus and to measure the systematic effects 
due to survey boundaries, variation of angular selec- 
tion function, galaxy biasing, and redshift space distor- 
tion. We adopt the halo-galaxy one-to-one correspon- 
dence (HGC) model of Kim et al. (2008) to assign mock 
galaxies to dark halos identified in a cosmological N-body 
simulation. For this purpose we made two N-body sim- 
ulations of the ACDM model having the WMAP 3 year 
cosmological parameters. The simulation and cosmolog- 
ical parameters adopted in the simulations (SI and S2) 
are listed in Table [5J 

The simulations were started at the initial redshift of 
Zi, and evolved to the present epoch after making Nstep 
global timesteps. Then dark matter particles at the 
present epoch are used to identify the dark halos. The 
CDM particles in high density regions are first found by 
using the Friend-of-Friend (FoF) algorithm with connec- 
tion length of 0.2 times the mean particle separation, 
and then gravitationally bound and tidally stable dark 
halos are identified within each FoF particle group (Kim 
& Park 2006). These dark halos include isolated, cen- 
tral and satellite halos composing FoF halos. The halos 
are required to contain at least 30 particles, which cor- 
responds to 2.5 X 10"/i"^Mq for 1024/1"! Mpc box size 
simulation (SI) and 6.8 x lO^/i-^Afo for lAii.&h-^ Mpc 
box size one (S2). The resulting halo samples include 
12,326,725 and 14,244,305 dark matter halos, and the 
corresponding mean halo separations are 4:Ah~^ Mpc and 
5.9/i~! Mpc, respectively. 

To select mock galaxies corresponding to those in 
BEST from the SI simulation, we sort the dark halos 
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Fig. 5. — A genus curve measured from the BEST sample (solid line with error bars), and its best-fit Gaussian eurve (long-dashed line). 
Curves writh smaller amplitudes were measured from the BEST sample of SDSS DR3 by Park et al. (2005). The observed genus curves in 
this plot have not been corrected for systematic biases. 



in mass and choose those above a mass cut so that the 
resulting halo set has mean halo separation of 6.1/i~^Mpc 
which is the mean separation between galaxies in BEST. 
These halos are identified as the galaxies that can be 
compared with those in BEST. 

We make 27 mock BEST samples originating at ran- 
dom locations within the SI simulation using exactly 
the same survey mask, angular selection function, ra- 
dial boundaries, and smoothing length, and also taking 
into account the peculiar velocities of the mock galax- 
ies by adding the line of sight component of the peculiar 
velocity to the distance of each galaxy, and we calculate 
the genus curve for each sample over a set of smoothing 
lengths. We do not enforce non-overlapping volumes for 
these mock samples, but the ratio of the simulation vol- 
ume to the BEST sample volume is approximately 40:1, 



so the 27 mock catalogs are approximately independent. 
The error bars in Figure [5] are the standard deviation 
of the genus curves from these 27 mock surveys. The 
derived statistics such as the amplitude (g), shift (Ai^), 
cluster abundance {Ac), and void abundance {Ay) pa- 
rameters are then measured from the genus curves. 

The genus calculated from the observational data suf- 
fers from systematic effects due to redshift space dis- 
tortion and complicated survey boundaries. To esti- 
mate the effects we use the difference between the genus 
curves measured from the real-space galaxy distribution 
in the whole simulation cube and from the redshift-space 
galaxy distribution in the mock surveys. For exam- 
ple, we measure the genus density g^ '^ at a smooth- 
ing scale Rg using all halos with mass Mh > 7.9 x 
IO^/i-^Mq {d = 6.1/i"iMpc) in the simulation cube. 
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Fig. 6. — Genus-related statistics as a function of Gaussian 
smoothing length Rq measured from the BEST sample of galaxies 
with Mr < —20.19. In the lower panels, a Gaussian field predicts 
Af = 0, Ac = Ay = 1. Open circles are those before the system- 
atic bias corrections. 

Halo positions in real space are used for g^^'^^^. The 
mean genus density gz™°^^ is also calculated from the 
27 mock SDSS surveys where the halos with Mh > 
7.9 X 10^^h~^MQ are observed in the simulation, with 
their rcdshift space positions. The difference between the 
two contains the information on the systematic effects of 
survey boundaries and rcdshift-spacc distortion (RSD). 
The observed genus in real space can be estimated by 
gr°^'^ — gz°^^gr™^^/gz™°'^^ assuming that the correction 
is model-insensitive. The correction factor is obtained 
from the ratio of the two statistics in the case of 5, Ac, 
and Ay, and from the difference in the case of lS.v. 

The correction is expensive because we need to find the 
correction factors for all statistics whenever we change 
the volume-limited sample or the smoothing scale. But 
there is a big advantage for presenting the observational 
data corrected for the survey boundary effects and RSD 
effects. It enables one to compare the observational data 
with theoretical predictions based on analytic calcula- 
tions or the analysis of the simulation data without tak- 
ing into account the particular survey characteristics. An 
empirical study on the systematic effects on the genus is 
presented in Appendix \K\ 

4.3. Scale- dependent Topology Bias 

To study the scale dependence of topology without the 
effects of changing large-scale structure in the survey vol- 
ume, we fix the sample definition (i.e. angular and radial 
boundaries) and change the smoothing length only. Fig- 
ureinishows the genus-related statistics (filled circles with 
error bars) as a function of Gaussian smoothing length 
in the case of the BEST sample and the data arc given 
in Table [3] Error bars arc obtained from 27 mock BEST 
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Fig. 7. — Genus curves measured from the BEST sample at throe 
smoothing scales, Rq = 6 (black line), 8 (red line), and lOh^^ Mpc 
(blue line). Their amplitudes are scaled by (RQ/dh^^Mpc)'^. 
These curves are not corrected for the systematic biases. 



surveys made in the SI ACDM simulation using the halos 
having the mean separation of 6.1/i~^Mpc. Open circles 
are the statistics before the systematic bias corrections 
are made. 

Figure [S] demonstrates that the deviation from the 
Gaussian-field topology is scale-dependent. The sec- 
ond panel shows that genus curve shifts towards in- 
creasingly negative values of Ai/ ("meatball" topology) 
as Rq increases from 6 to 10/i^^Mpc. The fraction- 
ally largest and statistically strongest deviations from 
the Gaussian field prediction arise for the void abun- 
dance parameter Ay, which rises from 0.7 at 6h~^ Mpc 
to 0.8 at 10/i~^Mpc (compared to Ay = 1). We caU 
these phenomena the scale-dependence of galaxy clus- 
tering topology. The cluster abundance Ac also dif- 
fers from the Gaussian-field prediction Ac = 1, but 
in a way that is approximately constant with smooth- 
ing length over this range. Note that the figures plot 
gRa^ = {G ohs / Vsampie) Rg^ , which is proportional to the 
genus per unit smoothing volume. 

To visually relate these measures to genus curve 
shapes, we show the genus curves for Re = 6, 8, and 
lO/i-^Mpc, scaled by {Rc/Qh-^Mpcf in Figure[71 

To understand the bias in the topology of galaxy clus- 
tering wc need to know the underlying matter field. We 
assume the matter field of the universe is equal to that 
of our ACDM simulation with WMAP 3 year parame- 
ters. The genus is calculated for the CDM field (of the 
full 1024/1^^ Mpc simulation cube) as in Figures 1 and 
2 of Park, Kim & Gott (2005) and compared to that 
of the observed genus in Figure [S] where deviations of 
gg/gm,Ac,glAc^m., and Av,g/Av,m from I, or deviation 
of Afg — Ai/'m from indicate the galaxy biasing with re- 
spect to matter. We find significant bias in the topology 
of the galaxy distribution for all four statistics. When 
compared to the matter distribution, galaxy clustering 
shows more complicated structures (higher g), has more 
meat-ball shifted topology (more negative Aj/), fewer 
number of superclusters (smaller Ac), and fewer num- 
ber of voids (smaller Ay). In particular, the void abun- 
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TABLE 3 

Genus-Related Statistics of the Scale-Dependence Samples in Figure [6] 



statistics 






«G 








6.0 


7.0 


8.0 


9.0 


10.0 


G 


410.0(373.8) ± 17.9 


270.3(247.0) ± 13.5 


184.7(168.1) ± 10.6 


133.7(121.8) ± 9.3 


97.5(89.4) ± 7.3 


Ai. 


-0.023(-0.027) ± 0.019 


-0.027(-0.027) ± 0.024 


-0.047(-0.045) it 0.030 


-0.068(-0.063) ± 0.033 


-0.072(-0.066) ± 0.036 


Av 


0.70(0.71) ± 0.04 


0.73(0.74) ± 0.05 


0.77(0.78) ± 0.06 


0.76(0.77) ± 0.06 


0.81(0.81) ± 0.07 


Ar: 


0.83(0.81) ± 0.05 


0.84(0.83) ± 0.07 


0.85(0.84) ± 0.08 


0.84(0.83) ± 0.08 


0.81(0.80) ± 0.07 



Note. Rq is the smootliing length in units of h^ Mpc. G is tlic amplitude of tlic genus curve, Av is tlie shift parameter, and A(^ and A^^ are cluster and void abundance 
parameters, respectively. All these values are systematic bias-corrected, and the observed values before tire systematic bias corrections are given irr parentheses. Uncertainty 
limits are estimated from 27 mocJt surveys in rodshift space. 
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sample (Mr < —20.19) arc used to measure the galaxy topology. 
The CDM distribution in our ACDM simulation is adopted as the 
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dance shows a very strong bias which has a strong scale- 
dcpendcncc. Any successful galaxy formation model 
should explain this topology bias in the distribution of 
galaxies. 

James et al. (2007, 2009) studied the scale depen- 
dence of topology using the 2dFGRS samples. James 
et al. (2009) measured the genus statistic as a function 
of scale from 5 to 8h~^ Mpc, using analytic predictions 
for both the linear regime (Gaussian random field) and 
weakly non-linear regime (perturbation theory; Matsub- 
ara 1994; Matsubara & Suto 1996) to detect the effects 
of non-linear structure formation on the genus statistic, 
and found that both Ay and Ac tend to fall below unity 
independently of scale and that neither of the analytic 
prescriptions satisfactorily reproduces the measurement. 

Quantitative meausures of the shape parameters of the 
genus curve before and after the bias correction shown 
in Figure [S] help us understand the degree of the RSD 
effects on the genus. We find that the survey boundary 
effects are very small because the volume-to-surface ra- 
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Fig. 9. — Amplitude of the genus curve in redshift space relative 
to that in real space from 27 SDSS mock surveys. Dashed line is 
the linear theory prediction of Matsubara (1996). 

tio is quite large for the SDSS DR7 data (see blue solid 
and black dashed lines in the top panel of Figure [18] in 
Appendix [Af ■ 

On 6/i^^Mpc smoothing scale the shape parameters 
{Ai^, Ac, and Ay) change only a little while the am- 
plitude reduction is as much as 8.8%. The RSD effects 
make the i^ > I part of the genus curve (Ac) have an 
amplitude lowered by 2.4%, and increase the v < —1 
part (Ay) by 1.4% relative to the best-fitting Gaussian 
genus curve (see Table 3). Even though the overall am- 
plitude of the genus curve is rather sensitive to the RSD 
effects, its shape is not. This fact is in agreement with 
the perturbative theoretical description of Matsubara & 
Suto (1996). 

To compare Matsubara (1996) 's linear prediction for 
the RSD effects with those calculated from simulation, 
we measure the amplitude of the genus curves of the 
27 mock survey samples both in redshift space and real 
space, which are denoted by G^(i/) and G^{v), respec- 
tively. Their ratios are plotted as a function of the 
smoothing length in Figure [9l Dashed line is the linear 
theory prediction according to Matsubara (1996) 



G'[v) _ 3V3 



Vu(l-u), 



where 



11 



i6/5)fb-' + (3/7)(/6^i)2 



(6) 



(7) 



3 1 + (2/3)/5-i + (l/5)(/fe-i)2- 
The b is the bias parameter and the function / is given 
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Fig. 10. — Genus-rclatcd statistics for subsets of two volume- 
limited luminosity samples with —19.0 > Afr > —21.0 and —19.5 > 
Mr > —21.5. Three connected points belong to the same volume- 
limited sample. The absolute magnitude ranges of the subsets are 
(-19.5, -20.1), (-19.8, -20.5), and (-20.1, -21.5) for the brighter 
sample, and (-19.0, -19.7), (-19.3, -20.2), and (-19.7, -21.0) for 
the fainter sample. These ranges are chosen so that each subset 
contains the same number of galaxies. All analysis use the same 
smoothing length of Rq = 5.9h~^ Mpc, and the systematic effects 
are corrected. Note that the error bars in this figure do not include 
cosmic variance — they are relevant for comparing the three points 
in each sample but not for comparing the samples to each other. 



by 

f{a,^, ^^) = d\nD/dhm « 0^„« + ^(1 + ^) (8) 

where D(t) is the linear growth factor and a is the 
scale factor. The linear theory predicts G^{v)/G'^{v) = 
0.991, while our rcsuh is 0.931 at Rg = 6/1"^ Mpc. 
C- {v) / G"^ {v) tends to slowly approach the linear pre- 
diction as the smoothing length increases. Matubara & 
Suto (1996) have pointed out that in weakly non-linear 
regimes the RSD tends to suppress the amplitude more 
strongly than the linear theory prediction, in agreement 
with our result. Figure 9 also shows that, apart from 
the sample-to-samplc fluctuation in the observed density 
field, the sample-to-sample variation in the large-scale 
peculiar velocity field causes about 1.5% uncertainty in 
the RSD correction on the Rq = 6/i~^ Mpc scale in the 
case of the SDSS DR7 sample. (Note that Fig. 9 shows 
the effects of RSD only with the density field in each 
sample fixed.) 

4.4. Dependence of Topology on Galaxy Properties 

Park et al. (2005) presented the first clear demon- 
stration of luminosity dependence of galaxy clustering 
topology: brighter galaxies show a stronger signal of 
"meatball" topology. To quantify the luminosity bias 
we use two volume-limited samples with —19.0 > M^ > 
—21.0 (corresponding to a comoving distance range of 
59.7 < R < 203.2/1-1 Mpc) and -19.5 > Mr > -21.5 
(59.7 < R < 253.0/1^1 Mpc). Each sample is divided 



into three subsets of galaxies belonging to three abso- 
lute magnitude ranges. To be free from confusion due to 
changing levels of discrete sampling of structures (ran- 
dom fluctuations in the number of galaxies that trace the 
density field), we choose the subsamples to contain the 
same number of galaxies. Since each of the three subsets 
includes the same number of galaxies distributed in the 
same volume of the universe, any difference in topology 
must be due to difference in luminosity. We then mea- 
sure the genus-related statistics to quantify the luminos- 
ity bias. The statistics measured for Rq = 5.9/i-iMpc 
arc plotted in Figure [TUl The scale is chosen because it is 
the mean separation of the galaxies in the subsets of the 
bright volume-limited sample. This smoothing length 
is also used for the subsets of the faint sample, whose 
d is 5.4/i-iMpc. The ^-coordinate of each filled circle 
is the median absolute magnitude of each subset, and 
three circles corresponding to three subsets of the same 
volume-limited sample are connected together. 

In Figure [10] there are two sets of connected points 
corresponding to the two volume-limited samples. Dif- 
ferences between the two sets are due to the fact that the 
samples enclose the universe with different outer bound- 
aries (i.e., cosmic variance) and should not be paid at- 
tention here. The luminosity dependence of topology 
appears within each volume-limited sample. 

To estimate the significance of the luminosity bias we 
make 20 mock subsets for each volume-limited sample 
by selecting 1/3 of the galaxies randomly without taking 
into account luminosity. Observational and mock sub- 
sets have the same sample volume, geometry, and galaxy 
number density. They differ only by galaxy luminosity. 
We use those mock subsets to estimate the error bars 
in Figure [TOl Therefore, they do not contain the cosmic 
variance, and show the significance of luminosity selec- 
tion only. 

It can be seen that the distribution of brighter galaxies 
tends to have a stronger "meatball" topology signature 
(more negative Aj^) and greater percolation of under- 
dense regions (lower Av). But the luminosity depen- 
dence for the shift (At/) parameter is somewhat mild, 
particularly for the faint sample. We do not detect any 
statistically significant dependence of g on luminosity. 

After luminosity we consider morphology. It is well- 
known that early-type galaxies are dominant in mas- 
sive clusters and late types are abundant in the field, 
and one naturally expects morphology dependence of 
the topology of galaxy clustering. The left panel of 
Figure [TT] shows the genus curves of early- (red curve 
with a smaller amplitude) and late-type galaxies with 
-20.99 <Mr < -20.19 at 9.1/i-iMpc Gaussian smooth- 
ing scale. We limit the absolute magnitude range to re- 
duce the luminosity-dependent bias and to inspect the 
morphology dependence of the genus only. 

We find the mean galaxy separations of the early- and 
late- type galaxy samples are 8.9 and 8.3 ft.-"'^Mpc, respec- 
tively. We discard the most edge-on late-type galaxies to 
make the mean galaxy separation in_the late-type sample 
equal to 8.9 /i^^Mpc. In Table IH d before the thinning 
out is given in parentheses. The smoothing length of 
Rg = 9.1/i-iMpc is chosen because it is the mean sep- 
aration of 'blue '-type galaxies, which we will use below 
to study the difference between morphology and color 
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Fig. 11. — Genus curves of early- and late-type galaxies (left panel) and red {u — r > 2.4) and blue (m — r < 2.4) galaxies (right panel) in 
the BEST sample with narrow luminosity range (—20.99 < Mr < —20.19) and at Gaussian smoothing scale of 9.1/i^^Mpc. These curves 
are not corrected for the systematic biases. 



dependence. The systematic effects are again corrected 
and error bars are estimated from 20 subsets drawn out 
of the volume-limited sample made by selecting galaxies 
randomly without considering rnorphology Each set has 
the mean galaxy separation of d = 8.9/i~^Mpc. It can be 
easily seen that the genus curve for early-type galaxies 
has a smaller overall amplitude, fewer voids, and more 
clusters, and is meat-ball shifted compared to that of 
late- type galaxies. These differences arc manifest in Fig- 
ure [1^ where wc show the genus-related parameters for 
early- and late-type samples separately. 

The upper panels are for the brighter sample men- 
tioned above, and the lower panels are for relatively 
fainter galaxies with -20.20 < M,. < -19.40. In the 
lower panels, the smoothing length Rq is 7.9 h~^Mpc, 
which corresponds to the mean separation of the early- 
type galaxies with magnitudes in this range (see Ta- 
ble [T2)) . Table |9] in Appendix [C] gives the data for the 
genus curves. Table |4] lists the measured genus-related 
statistics, both corrected and uncorrected (in parenthe- 
ses) for systematics. It is evident that significant mor- 
phology bias in the topology of the galaxy distribution 
exists, confirming the difference between the early- and 
latc-typc distributions predicted from hydrodynamic cos- 
mological simulations by Gott, Gen, & Ostriker (1996). 

Next, we examine whether the genus curve depends on 
galaxy color or not. Although the morphology and color 
are strongly correlated with each other, the roles of mor- 
phology and color in galaxy evolution are expected to 
be different. For example, galaxy color seems to depend 
mainly on the distance and morphology of the nearest 
neighbor galaxy, while morphology depends not only on 
neighbor but also on the distance to the nearest massive 
cluster of galaxies (Park & Hwang 2009). Thus, it is in- 
teresting to examine separately how the galaxy clustering 
topology depends on these two physical parameters. 

For the color comparison, the genus curves of red 
{u — r > 2.4) and blue {u — r < 2.4) galaxies with 
-20.99 < Mr < -20.19 are plotted in the right panel of 
Figure [TT] We create red and blue samples with the same 



number of galaxies by randomly throwing away some 
galaxies in the larger sample. Edge-on galaxies are dis- 
carded here. In Tablc[7]of Section 5.3, the mean separa- 
tion of galaxies in each sample and smoothing length are 
given. In this color-dependence analysis, we adopted the 
same sample used for the morphology-dependence study. 
The genus-related statistics for the color subsamples are 
shown in Figure [13] together with those of morphology 
samples for comparison. Table|9]in Appendix ICl gives the 
data for the genus curves, and Table[5]in Section 5.3 lists 
all the measured genus-related parameters. Similarly to 
the early-type galaxies, the distribution of red galaxies 
has a smaller overall amplitude of the genus curve, fewer 
voids, and more clusters, and is meat-ball shifted com- 
pared to that of blue galaxies. The result confirms the 
relative meatball shift between red and blue galaxies de- 
tected in a 2-D topology survey from the SDSS Early 
Data Release by Hoyle et al. (2002). 

There is one major difference between topology of mor- 
phology subsets and color subsets. The contrast in Ay 
between the early and late subsets is much larger than 
that between red and blue subsets. This can be also seen 
for fainter galaxies with —20.20 < Mr < —19.40 used 
in Section 5.3. The dependence of the void abundance 
on color is much reduced, while the dependence on mor- 
phology still exists in the fainter magnitudes as shown in 
the bottom right panel of Figure [121 

5. TEST OF GALAXY FORMATION MODELS 

5.1. Galaxy Formation Models 

Galaxies are the end product of non-linear gravi- 
tational evolution of primordial density fluctuations. 
Therefore, the spatial distribution of galaxies and its de- 
pendence on the internal properties of galaxies should 
depend both on initial conditions and on galaxy forma- 
tion and evolution processes. This enables us to test 
galaxy formation models as well as the models on pri- 
mordial fluctuations using topology. In this section we 
examine whether or not various models of galaxy forma- 
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TABLE 4 
Genus-related Statistics of the Morphology-Dependence Samples in Figure [121 





-20-19 > 


Mr 


> -20.99 


-19.40 > Mr 


> -20.20 




(59.7 < R < 


340 


.6h~^ Mpc) 


(59.7 < R < 242 


.2h~^ Mpc) 




Early 




Late 


Early 


Late 


Statistics 


d = S.9 




d ^ 8.9(8.3) 


d = 7-9 


d = 7.9(6.3) 


G 


128.1(117.4) ± 2.9 




141.7(129.8) ± 2.9 


66.6(61.2) ±2.1 


65.4(60.1) ±2.1 


Au 


-0.121(-0.111) ± 0.028 




-0.071(-0.061) ± 0.028 


-0.075(-0.067) ± 0.018 


0.012(0.020) ± 0.018 


^V 


0.68(0.67) ± 0.04 




0.85(0.84) ± 0.04 


0.83(0.85) ± 0.06 


0.96(0.99) ± 0.06 


^C 


0.93(0.91) ± 0.04 




0.76(0.74) ± 0.04 


1.09(1.03) ± 0.06 


0.84(0.80) ± 0.06 



ISfote 

Bystoi 



parentheses is the value before the thinning out the late- type 



Qplc. All genus-related statisti< 



tion are consistent with our measurement of topology of 
galaxy clustering, assuming that the difference between 
observation and models is entirely due to inaccuracy in 
the galaxy formation mechanism. This assumption is 
supported by the fact that the topology of the Luminous 
Red Galaxies in the SDSS analyzed by Gott et al. (2009) 
agrees very well with that of the mock galaxies identified 
in the ACDM model with the same cosmological param- 
eters (i.e. 3 year WMAP parameters) we adopt in this 
paper. Gott et al. measured the topology at very large 
scales of Rq = 21 and 34/i~^Mpc. In this almost linear 
regime, the results are sensitive only to the adopted cos- 
mological parameters. They use our simple HGC model 
to locate the mock LRGs in a large ACDM simulation. 
Gott et al.'s results indicate that our cosmological model, 
which adopts the initially Gaussian primordial fluctua- 
tions with the WMAP 3 year ACDM power spectrum, is 
consistent with the observation in the linear regime. If 
one finds a discrepancy between observation and models 
in the galaxy clustering topology on non-linear scales, 
therefore, it most likely arises from inaccuracy in the 
galaxy formation models. 



We adopt five galaxy allocation schemes applied to 
cosmological N-body simulations. The first is the HGC 
model described in Section 4.2. Each of the dark ha- 
los identified in our SI ACDM simulation of 2048^ CDM 
particles in a 1024/i~-'^Mpc size box is assumed to host 
one galaxy. The halos used in the HGC model can be 
either isolated or grouped, and in the latter case they 
can be the central one or satellites. The second is a 
Halo Occupation Distribution (HOD) model of Yang et 
al. (2007) refined by introducing a conditional luminos- 
ity function (van den Bosch et al. 2007) that accurately 
matches the SDSS luminosity function and the cluster- 
ing properties of SDSS galaxies as a function of their 
luminosity. The other three models that we examine 
are different implementations of semi-analytic models of 
galaxy formation (SAMs): Croton et al. (2006), Bower 
et al. (2006), and Bcrtonc et al. (2007). AU of them 
are set in the context of structure formation in a CDM 
universe as modelled by Millennium Simulation (Springel 
et al. 2005) of 2160^ particles in a 500/i"^Mpc size box 

( n^ = 0.25, r^A = 0.75, r^fc = 0.045, /i = 0.73, n^ = 1 

were assumed). However, they used different construe- 
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Fig. 13. — Genus-related statistics for the morphology {filled sym- 
bols) and color subsamples {open symbols) with a narrow luminosity 
range of —20.99 < Mr < —20.19. Open circles and squares arc for 
color subsets divided by the color cuts {u — r)cut = 2.4 and 2.6, 
respectively. 



tion of the dark matter merger trees and different im- 
plementation of the various physical processes involving 
the baryonic component of the universe. Croton et al. 
and Bower et al. invoked 'radio-mode' AGN feedback 
schemes to restrict gas cooling within relatively massive 
halos. The prescriptions in both models were set by the 
requirement that they approximately reproduce the ob- 
served break in the present day luminosity function at 
bright magnitudes, but their detailed implementations 
differ. Bertone et al. have used a variation of the SAM 
of de Lucia & Blaizot (2007), which is the later version 
of the galaxy formation models described in Croton et 
al., and implemented supernova feedback using a 'dy- 
namical' treatment of galactic winds (Bertone, Stoehr & 
White 2005). The treatment of fast recycling of ejected 
gas and metals predicts a relatively low abundance of 
dwarf galaxies in better agreement with observations but 
a rather larger number of bright galaxies than observed. 

5.2. Topology Test 

For a given luminosity cut of a volume-limited sample 
of the SDSS galaxies, the mean galaxy separation d and 
the maximum sample depth are determined. To compare 
the mock galaxies generated by the above five galaxy 
formation models in a fair way, we set the mean separa- 
tion between the mock galaxies in each model equal to 
the observed value by adjusting the lower limit of galaxy 
mass or luminosity. This gives us the relation between d 
and the limiting absolute magnitude Af^^cut- The outer 
boundary i?max of the volume-limited sample is deter- 
mined for the apparent magnitude limit of ?' = 17.6 and 
the absolute magnitude limit of Mr^cut- 



Comparisons among the observations and models are 
made for galaxies with the same number density or d. For 
example, we sort all galaxies (halos) in our HGC cata- 
log and find the mass cut above which the mean galaxy 
separation is equal to that of each volume-limited sam- 
ple of the SDSS galaxies. The resulting sample of mock 
galaxies with mass above the cut is to be compared with 
the corresponding observational sample. Mock galaxy 
subsets are similarly made for the HOD and three SAM 
galaxy samples. The mass or luminosity cuts are given 
in Table El 

The genus-related statistics measured from six volume- 
limited samples of the SDSS galaxies (dots with error 
bars) are plotted in Figure [14] at locations correspond- 
ing to the sample defining parameters (d, M,._cut, ^max) 
which are also given in Table [5] It should be noted 
that the observed volume-limited sample contains pro- 
gressively brighter galaxies and the sample depth i?max 
increases as the mean galaxy separation d increases. Fig- 
ure [T3] also shows the results for the mock galaxies of five 
galaxy formation models. Table [B] lists the statistics for 
all models and samples. 

Figure [15] shows the genus curves for the SDSS galaxies 
in the BEST sample and for the five sets of mock galaxies 
with 6.0 h~^ Mpc smoothing. Table [TOl in Appendix [Cl 
gives genus curves of the SDSS and mock galaxies. Four 
models, the exception being Bower et al. (2006), agree 
very well with one another for the genus density and shift 
parameters over the scales explored. Those four models 
agree well with the observations in terms of g, but Bower 
et al. (2006) agrees better for Ai/. The models give di- 
verse results for the cluster and void abundance parame- 
ters, though all of them predict Ay+Ac < 2 at all scales 
shown, indicating (Matsubara 1994) that effects beyond 
those of perturbative non-linear gravitational evolution 
must be involved. Bower et al. (2006) is again most dis- 
crepant with the observation in terms of Ac-, while HGC 
and HOD match the observed values almost perfectly. 
When Av is considered, HOD and all SAM models are 
inconsistent with the observations, while HGC is within 
the cosmic variance from the observation. Overall, no 
model reproduces all features of the observed topology. 
Only the Bower et al. (2006) SAM comes close to match- 
ing Av, but it fares worst with Ac and the genus density. 
Only HGC matches Ay, but it fails for Aj/. 

To estimate the statistical significance of the failure of 
each model we select two volume-limited samples with 
d = 6.10 and 7.llh^^ Mpc. We calculate the x^ statistic 



X 



E/,, model „,obs\2 /_2 



where the index i runs over the four genus-related statis- 
tics and j runs over the two volume-limited samples, erf.- 
is the variance of the i-th statistic calculated from 27 
mock surveys of the j-th volume-limited sample made 
in the SI simulation. A distribution of this x^ statis- 
tic is obtained from 64 mock volume-limited samples of 
the HGC galaxies made from the S2 simulation (see the 
histogram in Figure [16]) , where v°j'^ is replaced by the 
average value over the 64 mock samples. At each of the 
64 locations in the simulation two mock SDSS samples 
are made with depths of 340.6 and 388.3 h~^Mpc and 
with mass cuts of 10^^-^" and lO^^'^^/i^^M©, respectively. 
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Fig. 14. — Genus-related statistics riieasurcd from six volume-limited samples of SDSS galaxies (filled circles with error bars) compared 
with those of mock galaxies produced by five galaxy formation models. Resultsfor each model arc shown by a line connecting six symbols. 
The volume-limited samples are distinguished by the mean galaxy separation d, which also determines the absolute magnitude cut Afr.cut 
or the maximum sample depth i?max. Other sample definition parameters A4^r,cut and i?max arc given in the upper x-axes. 

TABLE 5 
Sample Defining Parameters of six volume-limited samples used in Figure [Til 



Rg 



^min ) ^n 



AfSDSS 
-""r.cut 



ArHGC 
^"cut 



A/HOD 



]\ fCroton 
-""r.cut 



i\ /■Bo\vcr 
-""r.cut 



5.2 


5.17 


59.7, 


288.1 


-19.80 


11.67 


-19.86 


-20.08 


-20.36 


-21.25 


6.0 


6.10 


59.7, 


340.6 


-20.19 


11.90 


-20.26 


-20.46 


-20.72 


-21.62 


7.2 


7.11 


59.7, 


388.3 


-20.50 


12.10 


-20.56 


-20.75 


-21.00 


-21.87 


8.0 


7.99 


59.7, 


422.2 


-20.70 


12.25 


-20.74 


-20.94 


-21.20 


-22.03 


9.2 


9.17 


59.7, 


458.6 


-20.90 


12.43 


-20.94 


-21.13 


-21.42 


-22.21 


9.9 


9.89 


59.7, 


477.8 


-21.00 


12.52 


-21.04 


-21.23 


-21.53 


-22.30 



Note. — Cols.(l) Gaussian smoothing length in units of h~^ Mpc; (2) Galaxy mean separation in units of h~^ Mpc; (3) Inner and outer 
boundary of the SDSS volume-limited samples in units of h~^ Mpc; (4) Absolute r-band magnitude at z = 0.1 of the SDSS sample {K 
and evolution corrected); (5) Logarithmic mass cut of halos in HGC model sample in units of h~^MQ; (6) Absolute r-band magnitude at 
2 = 0.1 (K and evolution corrected) of the HOD galaxy sample; (7), (8) and (9) Absolute rest frame r-band magnitude cuts of three SAM 
galaxy samples. 



Smoothing lengths of Rg = 6.0 and 7.2h~^ Mpc are ap- 
plied to these two mock samples, respectively. 

The histogram of x^ computed from the mock catalogs 
in this way is well described by a x^ distribution with 8 
degrees of freedom, as expected if the errors in the four 
statistics at the two smoothing lengths are essentially un- 
correlated. (Compare the histogram and dotted curve in 
Figure [m and see further discussion in Appendix |B|) As- 
suming this distribution and applying a 1-tailed x^ test, 
the probability for the HGC model to be consistent with 
the observation is only 0.4%. The HOD and three SAM 
models arc absolutely ruled out by this test, confirm- 
ing that the discrepancies seen in Figure [T3] are of high 



statistical significance. The least inconsistent model is 
HGC, followed by the HOD model, while the three SAM 
models are the most inconsistent. These discrepancies 
confirm the findings of Gott et al. (2008), using the DR3 
topology measurement and several N-body and hydrody- 
namic simulations, that existing galaxy formation models 
do not reproduce the topology of the SDSS main galaxy 
sample near the smoothing scales we are exploring. It is 
interesting to note that the HGC performs best among 
the five galaxy formation models even though it is con- 
strained by only one observational parameter, the mean 
galaxy number density. 

On the other hand, Gott, Choi, & Park (2009) found 
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TABLE 6 
Genus-related Statistics for the samples used in Figure [Til 



Statistics 


Rg 


SDSS 


HGC 


HOD 


SAM'^''°'™ 


g^^^Bowor 


gA^/[Bcrtono 


G 


5.2 


343.5(312.7) ± 18.5 


352.0 


357.6 


350.9 


328.9 


349.9 




6.0 


410.0(373.8) ± 17.9 


404.6 


409.4 


405.8 


380.9 


409.8 




7.2 


397.2(363.9) ± 20.9 


371.0 


377.0 


377.8 


351.6 


377.9 




8.0 


377.9(348.9) ± 18.6 


363.4 


369.6 


371.3 


345.1 


375.1 




9.2 


340.8(317.9) ± 15.3 


319.4 


331.6 


333.4 


315.6 


336.2 




9.9 


331.7(308.8) ± 14.7 


297.0 


311.2 


313.5 


297.7 


316.5 


Au 


5.2 


-0.005(-0.011)± 0.020 


0.049 


0.038 


0.029 


0.006 


0.034 




6.0 


-0.023(-0.027)± 0.019 


0.034 


0.034 


0.028 


0.000 


0.032 




7.2 


-0.033(-0.030)± 0.018 


0.023 


0.022 


0.025 


-0.006 


0.022 




8.0 


-0.060(-0.057)± 0.021 


0.009 


0.016 


0.012 


-0.017 


0.015 




9.2 


-0.084(-0.075)± 0.024 


-0.005 


-0.006 


-0.008 


-0.038 


0.006 




9.9 


-0.074(-0.069)± 0.025 


-0.008 


0.000 


-0.019 


-0.042 


-0.012 


Av 


5.2 


0.71(0.73) ±0.05 


0.75 


0.81 


0.81 


0.84 


0.88 




6.0 


0.70(0.71) ±0.04 


0.74 


0.81 


0.82 


0.84 


0.89 




7.2 


0.73(0.75) ±0.04 


0.75 


0.84 


0.83 


0.85 


0.91 




8.0 


0.79(0.81) ±0.04 


0.74 


0.84 


0.84 


0.84 


0.89 




9.2 


0.82(0.84) ±0.05 


0.74 


0.84 


0.83 


0.83 


0.87 




9.9 


0.81(0.82) ±0.05 


0.74 


0.83 


0.82 


0.79 


0.85 


Ac 


5.2 


0.80(0.79) ±0.05 


0.83 


0.82 


0.88 


0.94 


0.86 




6.0 


0.83(0.81) ±0.05 


0.88 


0.87 


0.93 


0.99 


0.90 




7.2 


0.93(0.91) ±0.05 


0.92 


0.93 


0.97 


1.04 


0.97 




8.0 


0.96(0.95) ±0.06 


0.95 


0.95 


0.99 


1.05 


1.01 




9.2 


1.00(0.98) ±0.06 


0.99 


0.98 


1.03 


1.08 


1.04 




9.9 


0.98(0.97) ±0.07 


1.00 


0.99 


1.04 


1.10 


1.06 



Note. Systematic corrected parameters Genus-related statistics measured from six volume-limited SDSS samples and from mock galaxy 
samples of five galaxy formation models. The smoothing length Rq is in units of h~^Mpc. The observed values are corrected for the 
systematic biases, and those before the correction are given in parentheses. 
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Fig. 15. — Genus curves for the SDSS galaxies (black solid line 
with error bars) in the BEST sample and five sets of mock galaxies. 
The curve for the observed sample is corrected for the systematic 
biases. 



that the topology of the luminous red galaxies (LRGs) is 
successfully fit by the HGC model we have tested here, 
which means that for the most massive galaxies, the HGC 
model docs seem to work well. Thus the formation of 
the most massive LRGs seems to be a cleaner problem 
that can be modeled with CDM simulations and sub- 
halo finding routines like those described by Kim & Park 
(2006). 

There are several implications that can be obtained 
from this test on how to fix these galaxy formation pre- 



FlG. 16. — Distribution (histogram) of the x^ statistic for the 
mock galaxies in the HGC galaxy assignment scheme. The x^ 
statistic measures the sample-to-sample fluctuation of the four 
genus-related statistics at two smoothing scales Rq = 6.0 and 
7.2h~^ Mpc (see the text for the definition of x^)- To obtain the x^ 
distribution 64 mock SDSS surveys are made in the S2 simulation. 
The distribution is fit well by a x^ distribution function (smooth 
solid line). Triangles indicate the x^ statistics between the galaxy 
models and the observation. 

scriptions to make them reproduce the observed galaxy 
clustering topology. The SAM models and HOD model 
should be modified so that the mock galaxies do not sepa- 
rate voids — connecting low density regions would lower 
Ay- The SAM models should also produce more con- 
nected massive clusters (i.e. smaller Ac). All models 
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should produce large-scale structure percolating at lower 
density levels (i.e. more negative /\v). The problems 
with the SAM models might be alleviated if the thresh- 
old interval for galaxy formation were more narrow (i.e., 
fewer galaxies in low mass halos), making voids cleaner 
and superclusters denser. 

In principle, of course, the discrepancies in Figure [15] 
could be a sign of non-Gaussian initial conditions, or of 
a linear power spectrum very different from the one we 
have adopted based on WMAP3. However, given that 
different galaxy formation prescriptions produce differ- 
ences comparable in magnitude to the observed discrep- 
ancies, even though no one model matches all aspects 
of the observed topology, we are inclined to ascribe these 
discrepancies to imperfect galaxy formation physics. The 
dependence of topology on galaxy morphology and color 
also supports this interpretation, as the differences be- 
tween the topology of early- and late- type galaxy subsets 
are larger than the discrepancies between the full galaxy 
sample and the observations. 

5.3. Topology Test for Color Subsets 

The SAM models give not only luminosity but also 
color of mock galaxies, which makes it possible to com- 
pare the observation and SAM models using the sam- 
ples divided by color. We adjust the lower limit of 
galaxy luminosity to make the mean galaxy separation d 
equal to that of the SDSS galaxies in the volume-limited 
samples with -20.2 < M,, < -19.4 (Sample L2) and 
-20.99 < Mr < -20.19 (Sample LI). Then each sam- 
ple is divided into red and blue subsets using a color cut 
{u — r)cut that makes the mean galaxy separation of each 
color subset equal to the corresponding color subsets of 
the SDSS sample. Table [7] lists the parameters defining 
the color subsets for the observed and simulated samples. 

Figure [T7] shows the genus-related parameters of the 
SDSS galaxies (filled circles with error bars) and mock 
galaxies (other symbols). The upper panels are for the 
relatively brighter sample of galaxies in LI, and the 
lower ones are for those in L2. We show two sets of 
error bars for the observational data points. The shorter 
ones should be used when difference between blue and 
red color subsets is a matter of concern, but the longer 
ones, to which the cosmic variance is included, should be 
looked at when a model is compared with the observa- 
tion. 

It can be seen in Figure [T7] that SAM predictions of 
the color-dependence of topology are often qualitatively 
correct but incorrect in magnitude. Failure of the SAM 
models in reproducing the observed clustering topology 
of color subsets is more evident for the bright blue galax- 
ies of Croton et al. (2006) in the upper left panel of 
Figure 1171 They have too high genus density and their 
difference in g with red galaxies is too large. All three 
models predict relatively too positive Ai/ for both red 
and blue galaxies compared with the combined sample. 
The right panels of Figure [17] show that Ay of both red 
and blue mock galaxies is again too large in all three 
SAM model compared with the combined sample. In the 
case of Bertone et al. and Croton et al., the Ay and Ac 
parameters of the brighter sample (upper right panel) are 
nearly the same for red and blue subsets, which clearly 
contradicts the observations. 

To summarize, the SAM models fail to reproduce the 



observed clustering topology of galaxies divided accord- 
ing to color. Their problems for the combined sam- 
ple persist in the color subsets, and new problems are 
added. In particular, disagreement with the observation 
is more serious for bright galaxies. There should be more 
bright blue galaxies in superclusters so that clusters are 
more connected in the distribution of bright blue galax- 
ies. Both bright blue and red galaxies should be formed 
less frequently in void regions so that voids can be more 
connected with one another. The threshold for formation 
of both red and blue galaxies should be at lower density 
on average so that the percolation of large-scale structure 
occurs at lower density. 

6. SUMMARY 

We use the SDSS DR7 main galaxy catalog supple- 
mented with missing rcdshifts and with increased spec- 
troscopic completeness to measure the galaxy clustering 
topology over a range of smoothing scales. The distribu- 
tion of galaxies observed by the SDSS reveals extremely 
diverse structures. 

A volume-limited sample, BEST defined by Mr < 
—20.19, enables us to measure the genus curve with 
amplitude of G = 378 ± 18 at a smoothing scale of 
6h~^Mpc, with the quoted uncertainty including all sys- 
tematics and cosmic variance. The amplitude is 5.4 times 
larger than our previous measurement using the SDSS 
DR3 sample (Park et al. 2005) and the uncertainty 
decreases from 10.5% to 4.8% at the same smoothing 
length of Rg = 6/i~^Mpc. We calculate the galaxy clus- 
tering topology over the interval from Rq = 6/i~^Mpc 
to 10/i~^Mpc, and find mild scale-dependence for the 
shift (Al') and void abundance (Ay) parameters. The 
measured genus curve is qualitatively similar to the form 
predicted for Gaussian primordial fiuctuations (Hamilton 
et al. 1986), but the differences are statistically signifi- 
cant at these scales: a shift of the peak towards negative 
ly, and fewer isolated voids and isolated clusters than 
the Gaussian prediction {Ac and Ay < 1). The bias 
in topology of galaxy clustering with respect to that of 
matter is measured by assuming that the matter density 
field is given by our ACDM N-body simulation. We de- 
tect strong topology bias in galaxy clustering, which is 
also scale-dependent. 

We confirm the luminosity dependence of galaxy clus- 
tering topology discovered by Park et al. (2005). The 
distribution of brighter galaxies is more shifted towards 
"meatbair topology (lower Af) and shows greater per- 
colation of voids (lower Ay). We find galaxy clustering 
topology depends also on morphology and color. Even 
though early (late)-type galaxies show topology similar 
to that of red (blue) galaxies, morphology-dependence of 
topology is not identical with color-dependence. In par- 
ticular, the void abundance parameter Ay depends on 
morphology more strongly than color. 

We tested five galaxy formation models, which are 
used to assign galaxies to the outputs of N-body sim- 
ulations. Three of these are semi-analytic models, one is 
an HOD model that assigns galaxies to halos with pa- 
rameters tuned to match other clustering statistics, and 
one is a scheme that assigns galaxies to halos and sub- 
halos. None of them reproduces all aspects of the ob- 
served topology, though the differences from one model 
to another are comparable to the discrepancies with the 
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TABLE 7 
Sample defining parameters of the color-dependence samples in Figure [Tt] 



Sample LI 



Sample L2 



Samples 



Luminosity 



(u ~ r)cut 



i?G 



Luminosity 



ju - r)cut d Rg 



SDSS 
Croton et al. 
Bower et al. 
Bertone et al. 



-20.19 > Mr > -20.99 
-20.46 > Mr > - 21.21 
-20.71 > Mr > - 21.51 
-21.60 > Mr > - 22.28 



2.40 


9.1 


9.1 


-19.40 > Mr > 


-20.20 


2.40 


6.9 


7.0 


1.97 


9.1 


9.1 


-19.67 > Mr> 


-20.44 


2.18 


6.9 


7.0 


1.15 


9.1 


9.1 


-20.00 >Mr> 


-20.68 


1.20 


6.9 


7.0 


1.30 


9.1 


9.1 


-20.81 > Mr > 


-21.59 


1.77 


6.9 


7.0 



Note. Cols. (2) and (5) Absolute r-band magnitude limit at z = 0.1 for the SDSS sample (K and evolution eorrected) and Absolute rest 
frame r-band magnitude limit for the SAM models; (3) and (7) u — r color cut used for the color subsets; (4) and (8) mean separation in 
units of h~^ Mpc between galaxies in the sample; (5) and (9) smoothing length in units of h~^ Mpc. 



TABLE 8 

Genus-related statistics for the color- dependence samples in Figure [TtI 







Sampl 


eLl 








Sample 


L2 


Statistics 


Red 






Blue 




Red 




Blue 


SDSS samples 


G 


126.6(116.2) 


±2.9 


138.1 


3(127.3) : 


±2.9 


95.6(86.2) ±2.7 


101.4(91.4) ±2.7 


Au 


-0.091(-0.091) 


±0.016 


-0.047(-0.047) 


±0.016 


-0.056(-0.042) 


±0.016 


0.025(0.039) ± 0.016 


Av 


0.73(0.73) ± 0.04 


0.81 


(0.82) ±0.04 


0.82(0.86) ±0.05 


0.83(0.87) ±0.05 


Ac 


0.98(0.93) ±0.04 


0.79(0.76) ± 0.04 


1.01(0.97) ±0.04 


0.80(0.77) ±0.04 


Croton et al. 


G 


161.0 






124.7 




99.9 




80.6 


Au 


0.041 






0.002 




0.061 




-0.050 


Av 


0.84 






0.87 




0.93 




0.85 


Ac 


1.06 






1.04 




0.87 




1.15 






Bower 


et al. 












G 


140.5 






120.7 




92.4 




81.2 


Av 


0.004 






-0.027 




0.044 




-0.032 


Av 


0.94 






0.85 




0.92 




0.85 


Ac 


0.96 






1.12 




0.84 




1.16 


Bertone et al. 


G 


148.8 






127.6 




92.4 




84.4 


Av 


0.024 






0.021 




0.023 




-0.019 


Av 


0.94 






0.92 




0.99 




0.87 


Ac 


1.02 






1.01 




0.82 




1.07 



Note. Genus-related statistics for the color subsets of the SDSS samples and the mock samples of three SAM models. The observed values 
before the systematic bias corrections are in parentheses. 



observations. For this reason, and because the initially 
Gaussian ACDM model successfully reproduce the ob- 
served topology of LRGs at large scales, we attribute the 
discrepancies to failures of the galaxy formation model 
rather than non-Gaussian initial conditions. The semi- 
analytic models can also predict the topology of color 
subsets, but none of them fully captures the observed 
topology differences between red and blue galaxies. 

In future work, we will investigate models with non- 
Gaussian initial conditions to see what levels of primor- 
dial non-Gaussianity can be ruled out by our measure- 
ments. In principle, the high-precision topology mea- 
surements presented here and by Gott et al. (2009) can 
provide valuable constraints on non-standard inflation- 
ary models or alternative hypotheses for the origin of 
primordial fluctuations. 

Appendix A details our estimates of systematic biases 
in the genus curve measurements, demonstrating that the 
dominant effect is peculiar velocity distortions in redshift 
space rather than sample geometry or boundary effects. 
Except where noted otherwise, observational measure- 
ments in this paper are corrected for these biases, so they 



can be compared to theoretical predictions in real space 
with periodic boundaries. Appendix B investigates error 
covariances, showing that while the individual points on 
the genus curve have strongly covariant errors, the statis- 
tics G, Aj/, Av, and Ac are approximately independent. 
Appendix C tabulates the full genus curves for our best 
samples, complementing the statistics recorded in earlier 
tables. 

While future surveys will use luminous galaxies and 
emission-line galaxies to probe structure in the distant 
universe, the SDSS DR7 sample is likely to remain the 
definitive map of large scale structure at low redshift 
traced by a broad spectrum of galaxy types, for the fore- 
seeable future. The measurements in this paper char- 
acterize the topology of this definitive sample, attaining 
unprecedented statistical precision and providing a valu- 
able test for future models of primordial fluctuations and 
galaxy formation physics. 

We thank Xiaohu Yang for providing us with the HOD 
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APPENDIX 
EFFECTS OF SURVEY SYSTEMATICS ON THE GENUS GURVE 

In this paper we present the observed genus-related statistics corrected for the systematic effects. We aim to obtain 
the statistics in real space with no boundary and angular selection effects. Such results can be readily compared 
with theoretical predictions calculated from mock galaxies simulated in a variety of cosmological models without going 
through a full analysis taking into account the peculiar velocities and the complicated angular and radial survey 
selection functions of a particular survey. To understand the effects of the systematics on the genus curve step-by- 
step, we first calculate the genus curve from the number density field of the halos with Mh > 7.9 x 10^^/i~^Mq 
{d = 6.1/i~^Mpc) in the whole simulation cube (SI simulation) using periodic boundary conditions. The halos are 
the ones identified in our ACDM simulation, which consist of the isolated halos, the central halos, and the satellite 
subhalos. Each of the halos is assumed to contain one galaxy. 

The blue solid line in the top panel of Figure [TSl is the corresponding genus curve when Rq = 6/i~^Mpc. This is 
the true genus curve we hope to measure with no RSD effects and no survey boundary effects (but with some shot 
noise effects). We then make 27 mock surveys of these 'galaxies' in the simulation cube with radial boundaries and 
angular selection function identical to SDSS DR7 trimmed as shown in Figure 1. The mock galaxies are located at 
their real-space positions, and the resulting mock survey samples are free of the RSD effects. The short-dashed line in 
the top panel of Figure [TH] is the genus curve averaged over these 27 samples. It can be seen that the radial boundary 
and angular selection effects introduce little change in the genus curve in the case of the main part of SDSS DR7 
thanks to the large volume-to-surfacc ratio and high and uniform angular selection function of the sample. A large 
bias in the genus curve is generated when galaxies are observed in redshift space. The red line in the top panel of 
Figure [18] is the case when both RSD effects and the effects due to survey boundaries and angular selection function 
are taken into account in the 27 mock SDSS DR7 surveys. The difference between red line and dashed line is due to the 
RSD effects. The amplitude of the genus curve is decreased mostly due to the smoothing of structure by small-scale 
peculiar velocities of galaxies, but its shape does not change much. The genus curve near the median- volume threshold 
{ly = 0) is slightly pushed to the negative threshold direction, thus making Ai/ decrease. Both the number of clusters 
and voids are decreased. The number of high density regions is decreased more because massive clusters of galaxies 
are strongly clustered with one another and two clusters often appear connected along the line-of-sight in redshift 
space. On the other hand, a void can be affected when strong fingers-of-God protrude into the void from surrounding 
clusters in redshift space. The number of voids is decreased only a little between the blue and red lines in the top 
panel of Figure [TSl because it is difficult for voids identified on Rq = 6ft.~^Mpc scale to be erased by this process. But 
the number of voids relative to those expected from the best-fitting Gaussian curve (Ay) is rather increased. 

It is also interesting to know how the genus curve depends on the type of the density tracer. A comprehensive 
study on this issue has been made by Park, Kim & Gott (2005), who used a set of N-body simulations of AGDM and 
the Standard Gold Dark Matter Universe to examine the dependence of the genus curve on gravitational evolution, 
biasing, RSD, smoothing scale, and cosmology. In the middle panel of Figure [TH] we show how the genus curve changes 
at Rg = 6/i~^Mpc scale as the density tracer changes from dark matter particles to dark halos with different mass 
limits. The genus curve for the halos with Mh > 7.9 x 10^^/i~^Mq {d = 6.1/i~^Mpc) is plotted by a blue solid line 
again. The green solid line is the genus curve for the less massive halos with mass Mh > 2.5 x 10^^ h~^ Mq that are 
sparsely sampled to match the mean galaxy separation to 6.1/i~^Mpc. The green dashed line is for the full less massive 
halo data with d = 4.4/i^^Mpc. The figure demonstrates the dark matter density field (black solid line) and the halo 
number density fields (blue solid line) have very different topology at 6ft.~^Mpc Gaussian smoothing scale. The Ai^ 
parameter is positive for both matter particles and halo. The Ay and Ac parameters for matter are close to 1 and 
much larger than those of halos. For the Ay parameter, the difference is much larger. The strong bias in the observed 
clustering topology of the SDSS galaxies with respect to that of matter is shown in Figure [H On the other hand, the 
relatively more massive halos with Mh > 7.9 x 10^^h~^MQ describe fewer superclusters or fewer voids at fixed volume 
fractions (smaller Ay and Ac) when compared with that of halos with Mh > 2.5 x lO^^/i~^M0. 

COVARIANCE MATRIX OF THE GENUS-RELATED STATISTICS 

To understand how the measurements of the genus at different levels are correlated with one another, in Figure [12] 
we show the covariance matrix 






(Bl) 



where SGi = G{i'i) — Giyi) is the difference of the genus from the mean at a threshold z/^, and ai is the standard 
deviation of the genus at Vi. Estimation is made using two AGDM simulations, SI and S2. (See Table[5]for simulation 
parameters. Detailed explanation of these simulations is given in section 4.2.) 

The bottom panel of Figure [TO] shows the covariance measured from the halos with d = 6.1/i~^Mpc. The whole 
simulation cube of the S2 simulation with size of 1433.6/i~^Mpc is divided into 64 subcubes of size 358.4/i~^ Mpc, and 
the 64 genus curves from these subcubes arc used to calculate the matrix. The filled and open circles represent positive 
and negative covariances, respectively. The covariance along the diagonal line is one by definition, and the radius of 
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Fig. 18. — Tests for systematic effects on the genus curve obtained from the HGC galaxies in the SI simulation, (top panel) Effects of 
the SDSS survey mask and the redshift space distortion. The blue solid line is obtained using all HGC galaxies in the whole simulation 
cube at their real space positions. The dashed line and thin solid line are the genus curve averaged over 27 mock SDSS surveys of the HGC 
galaxies in real and redshift spaces, respectively, (middle panel) Effects of different tracers on the genus curve. The blue solid line is again 
the real space genus curve from all HGC galaxies with Mj, > 7.9 X 10^^h~^MQ and d = 6.1/i~^Mpc in the whole simulation cube. The 
green solid line is for those with M^ > 2.5 X 10^^h~^MQ sparsely sampled to have d = 6.1/i~^Mpc. The genus curve of the galaxies with 
Mfi > 2.5 X lQ^^h~^MQ and with the full sampling {d = 4.4h~^Mpc) is shown by a green dashed line. The difference shows the effects of 
shot noise. The black solid line is the genus curve of the dark matter particle distribution, (bottom panel) The genus curves of the SDSS 
DR7 BEST sample before (thin red line) and after (blue line with error bars) systematic bias corrections. A Gaussian smoothing length of 
Rq = 6.0/i~^Mpc is used in all cases. 



circlesare proportional to the covariance. The covariance plot in the upper panel is obtained from the matter particles 
with d = 0.5/i^^ Mpc of the SI simulation divided into 512 subcubes of size 128.0/i~^ Mpc. In both plots a smoothing 
length of Rg = 6.0/i"^ Mpc is used. 

The covariance seems to be smallest for threshold levels z^ = 0, and ±-\/3, where the extrema of the genus curve 
are located (see Gott et al. 1990 for a similar finding for i/ = level). There is a positive covariance strip along the 
diagonal with a width of about v3/2 with the width of positive and negative correlation ranges slightly depending on 
V. There are also two negative covariance strips shifted from the diagonal by ±v3, two positive ones with shifts of 
±2-^3, and so on. The magnitude of the covariance decreases rather slowly from the diagonal. 

To understand how the genus-related statistics employed in this paper are correlated with one another, we calculate 
the covariance matrix 



< Au.Am, > 



SiSj 
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Fig. 19. — Covariance matrix of the genus between two threshold levels Ui and Uj. (upper panel) The SI simulation cube is divided 
into 512 subcubes of size 128.0/i~^ Mpc, and the number density field of dark matter particles in the subcubes are used to obtain 512 
genus curves and the covariance matrix using Rq = 6.Qh~^ Mpc. (lower panel) The S2 simulation cube is divided in 64 subcubes of size 
358. 4h~^ Mpc, and the covariance matrix is obtained from the 64 genus curves of the halo distribution in the subcubes. Smoothing length 
of Rg = 6.0h~^ Mpc is used. Covariance is proportional to the size of circle, and open circles indicate negative covariances and filled circles 
indicate positive covariance. 

where Ui is the i-th one of the four genus-related statistics, Am = Ui — Ui^ and Ui is the mean of the i-th statistic 
calculated from 512 genus curves of the dark matter particles in the 512 subcubes of size 128/i~^Mpc taken from the SI 
simulation. The smoothing length used is Rg = 3ft.~^Mpc. The average in Equation (B2) is taken over 512 subcubes, 
and Si is the standard deviation of the i-th statistic from its mean. The matrix is symmetric, and its diagonal elements 
are one by definition. Figure 19 shows that the genus-related statistics are roughly independent of one another, with 
the magnitude of the covariance ranging from 0.08 to 0.31. Thus, it is roughly acceptable to test a model using a 
simple x^-test with these four statistics and the quoted observational errors, ignoring covariance, a conclusion further 
supported by the mock catalog tests shown in Figure [121 

TABLE OF GENUS CURVES 

In this section, we give the tables of the genus curves for the convenience of readers who wish to use galaxy clustering 
topology to test their galaxy formation or cosmological models. In Table [5] the genus is given as a function of volume- 
fraction threshold level for the observational morphology and color subsamplcs. The genus curves for the SDSS galaxies 
in the BEST sample and the five sets of mock galaxies measured in Section 5.2 are given in Table ITOl The genus values 
of the SDSS galaxies are corrected for the systematic biases, and the uncorrected values are in parentheses. These two 
curves are plotted in the bottom panel of Figure [TH Electronic forms of these tables are available from the authors 
on request. Models can also be tested against the genus curve statistics reported in earlier tables. 
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TABLE 9 

The Genus values at a given threshold level of the Morphology and Color Samples in Figure [T2I and [TtI 
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-7-0 


-13.7 


-12.6 


-9-7 


-7.1 


-2.9 


-17.8 


-17.1 


-9-4 


-7-2 


-20.5 


-17.3 


-16.5 


-10.0 


-2.8 


-17.4 


-19.6 


-12.2 


-14.3 


-24.2 


-21.8 


-18.8 


-12.4 


-2.7 


-21.0 


-24.3 


-14-8 


-15.5 


-23.9 


-28.6 


-15.7 


-13.7 


-2.6 


-22.8 


-28.5 


-14.3 


-17.8 


-28.1 


-34.8 


-18.4 


-18.4 


-2.5 


-30.4 


-32.8 


-14.6 


-19.8 


-30.2 


-35.8 


-15.8 


-20.2 


-2.4 


-29.2 


-36.6 


-18-7 


-20.9 


-32.5 


-38.9 


-22.6 


-23.6 


-2.3 


-35.4 


-43.0 


-21-8 


-21.9 


-44.1 


-39.6 


-25.1 


-26.4 


-2.2 


-38. 8 


-46.4 


-21-0 


-21.9 


-43.1 


-42.8 


-28.4 


-26.2 


-2.1 


-42.7 


-47.2 


-19-0 


-26.3 


-40.3 


-42.4 


-32.7 


-36.6 


-2.0 


-36.7 


-47.4 


-22.6 


-25.9 


-38.7 


-43.2 


-37.4 


-35.8 


-1.9 


-40.1 


-40.1 


-18.8 


-25.0 


-34.4 


-48.3 


-38.4 


-38.0 


-1.8 


-39.2 


-43-1 


-26-1 


-32.4 


-37.5 


-50.4 


-34.8 


-42.1 


-1.7 


-37.1 


-54.6 


-25-1 


-29.1 


-40.5 


-45.7 


-41.4 


-44.5 


-1.6 


-39.4 


-57.9 


-23-5 


-20.2 


-37.7 


-40.9 


-35.9 


-31.1 


-1.5 


-26.5 


-48.2 


-24.7 


-21.5 


-37.4 


-39.7 


-29.1 


-27.1 


-1.4 


-23.9 


-44.0 


-24.8 


-23.1 


-36.3 


-36.9 


-23.5 


-29.3 


-1.3 


-20.8 


-34.9 


-20.4 


-24.6 


-22.6 


-43.1 


-14.0 


-22.9 


-1.2 


-16.5 


-34.4 


-4.4 


-11.7 


-11.0 


-35.6 


-13.2 


-19.4 


-1.1 


3.2 


-27.7 


4.3 


-4.1 


4.7 


-21.2 


-19.7 


-10.4 


-1.0 


24.3 


-5-2 


5-7 


-1.5 


32.1 


2.1 


1.5 


-0.3 


-0.9 


55.0 


17.9 


25-9 


1.0 


28.4 


26.3 


31.8 


18.6 


-0.8 


65.0 


49.7 


28.3 


16.5 


46.7 


31.6 


43.1 


21.1 


-0.7 


79.5 


57-3 


41-2 


30.4 


78.7 


45.9 


43.6 


32.6 


-0.6 


100.7 


83-1 


49-6 


30.6 


91.1 


69.9 


53.5 


47.3 


-0.5 


112.1 


92-9 


37-8 


37.3 


105.7 


78.6 


77.4 


61.6 


-0.4 


123.0 


82-0 


48-3 


36.8 


124.3 


108.9 


79.9 


68.1 


-0.3 


123.6 


121-5 


62-0 


41.9 


124.0 


137.4 


85.5 


86.4 


-0.2 


113.3 


151-4 


68-9 


52.8 


124.7 


146.1 


90.4 


100.6 


-0.1 


118.9 


173-3 


63-9 


65.5 


111.2 


153.0 


94.8 


98.5 


0.0 


134.5 


161-1 


68-6 


65.2 


127.0 


150.1 


97.6 


102.8 


0.1 


119.0 


157.0 


59-2 


68.9 


111.5 


143.5 


95.5 


110.6 


0.2 


111.7 


118-4 


58-1 


66.8 


114.3 


128.1 


74.5 


96.6 


0.3 


85.5 


112.1 


54-6 


68.9 


81.4 


100.9 


69.6 


83.9 


0.4 


77.7 


96.4 


51-0 


65.2 


81.9 


104.2 


78.4 


84.9 


0.5 


73.8 


84.3 


40-1 


56.7 


72.5 


85.0 


62.3 


76.9 


0.6 


52.9 


60.8 


33.7 


35.4 


60.1 


55.0 


48.8 


63.1 


0.7 


44.2 


32.8 


25.7 


29.3 


49.1 


44.2 


33.4 


34.9 


0.8 


22.1 


7-5 


9-5 


7.9 


28.1 


15.5 


21.1 


28.9 


0.9 


6.9 


-lS-7 


4-0 


-6.8 


1.9 


-11.9 


12.1 


14.6 


1.0 


-12.7 


-32-6 


0-1 


-10.7 


-20.6 


-20.2 


-4.2 


-1.0 


1.1 


-32.1 


-35.0 


-14.2 


-10.0 


-35.5 


-22.1 


126.7 


90.8 


1.2 


-44.2 


-48.6 


-31.4 


-15.4 


-37.6 


-30.4 


-36.4 


-22.8 


1.3 


-57.6 


-36-1 


-37-1 


-26.5 


-49.3 


-26.8 


-40.0 


-20.6 


1.4 


-50.6 


-30-9 


-31-1 


-25.1 


-45.5 


-37.9 


-45.4 


-31.9 


1.5 


-47.4 


-34-5 


-30-3 


-23.8 


-53.5 


-45.7 


-44.0 


-30.6 


1.6 


-46.0 


-37-7 


-33-5 


-25.1 


-70.6 


-53.9 


-40.3 


-35.7 


1.7 


-50.7 


-44.9 


-30-2 


-26.7 


-56.6 


-43.3 


-38.1 


-38.5 


1.8 


-50.5 


-46.1 


-23.6 


-23.5 


-45.1 


-42.7 


-29.4 


-39.7 


1.9 


-40.5 


-52.7 


-25.9 


-20.2 


-38.0 


-43.5 


-33.5 


-27.3 


2.0 


-35.9 


-46.2 


-25.7 


-16.4 


-36.6 


-48.6 


-32.0 


-29.6 


2.1 


-37.0 


-42.0 


-20-1 


-16.1 


-42.7 


-42.9 


-31.3 


-30.4 


2.2 


-38.1 


-32.5 


-19-7 


-14.9 


-34.3 


-40.0 


-32.5 


-31.3 


2.3 


-33.6 


-35.8 


-14.9 


-14.3 


-33.9 


-35.6 


-25.6 


-21.6 


2.4 


-28.7 


-30.7 


-15.2 


-12.1 


-29.7 


-29.2 


-20.0 


-17.6 


2.5 


-24.8 


-26.1 


-11.5 


-11.7 


-25.9 


-26.4 


-16.0 


-15.7 


2.6 


-19.4 


-22.8 


-13.6 


-10.0 


-18.2 


-24.6 


-16.9 


-15.6 


2.7 


-17.5 


-20.3 


-11-7 


-9.0 


-14.1 


-21.7 


-15.7 


-16.2 


2.8 


-12.0 


-18.4 


-9-8 


-6.9 


-13.3 


-16.6 


-11.8 


-13.1 


2.9 


-8.6 


-11.9 


-S-8 


-4.3 


-11. S 


-11.5 


-7.4 


-9.8 


3.0 


-7-3 


-11-0 


-5-6 


-4.6 


-8.2 


-10.7 


-5.9 


-3.1 



Note. All the values are systematic bias-corrected. 
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TABLE 10 

The Genus at a given threshold level of the Samples in Figure [Tsl 



u 


G( 


'')SDSS 


*^('^)hGC G('^)Bcrtoiic 


'^(^)hod 


<3(^)Croto„ <3(>')Bowor 


-3.0 


-37.2(-33. 


3) ± 6.6 


-27.0 


-31.0 


-30.2 


-30.4 


-36.8 


-2.9 


-43.1{-37. 


3) ± 7.3 


-35.4 


-38.0 


-37.0 


-35.5 


-43.0 


-2.8 


-53.3(-45. 


5) ± 7.9 


-43.1 


-45.4 


-43.8 


-46.9 


-51.0 


-2.7 


-67.9(-62. 


5) it 8.2 


-50.9 


-54.3 


-52.8 


-55.1 


-60.2 


-2.6 


-72.2{-65. 


3) ± 8.8 


-60.0 


-62.3 


-61.9 


-65.6 


-70.1 


-2.5 


-83.2(-76.1 


± 10.8 


-70.8 


-70.5 


-72.2 


-77.9 


-80.2 


-2.4 


-96.5(-SS.8 


± 10.8 


-82.6 


-84.7 


-82.9 


-90.9 


-90.9 


-2.3 


-93.9(-SS.l 


± 10.1 


-93.4 


-97.1 


-95.8 


-98.9 


-105.7 


-2.2 


-107.0(-99.5 


± 11.3 


-105.4 


-116.2 


-109.8 


-109.6 


-118.7 


-2.1 


-104.5(-98. 


2) ± 9.8 


-116.1 


-129.2 


-125.7 


-125.9 


-133.5 


-2.0 


-118.6(-112.7 


± 10.8 


-124.1 


-145.2 


-141.1 


-136.6 


-146.5 


-1.9 


-135.0(-124.0 


± 10.5 


-131.6 


-149.5 


-149.5 


-149.7 


-157.8 


-1.8 


-132.2(-122.1 


± 12.1 


-136.0 


-152.0 


-153.9 


-151.0 


-168.1 


-1.7 


-130.1(-123.7 


± 11.6 


-136.6 


-148.1 


-151.8 


-148.1 


-163.9 


-1.6 


-121.6(-115.2 


± 15.5 


-131.0 


-145.4 


-145.0 


-136.4 


-161.7 


-1.5 


-116.3(-10S.l 


± 12.5 


-120.0 


-134.3 


-132.1 


-131.0 


-153.4 


-1.4 


-90.8(-84.3 


± 11.7 


-104.4 


-115.8 


-112.1 


-110.0 


-129.4 


-1.3 


-78.9(-72.2 


± 11.9 


-83.3 


-89.7 


-96.5 


-84.3 


-104.7 


-1.2 


-47.7(-43.1 


± 10.3 


-55.5 


-60.0 


-69.3 


-52.2 


-74.0 


-1.1 


8.0(7.6 


± 10.6 


-24.5 


-31.4 


-31.0 


-17.2 


-37.6 


-1.0 


19.7(19.3 


± 12.1 


14.8 


6.1 


10.0 


13.3 


6.9 


-0.9 


66.7(57.2 


± 12.6 


59.4 


60.0 


58.6 


63.3 


48.7 


-0.8 


112.0(101.1 


± 15.7 


102.9 


112.1 


105.1 


109.2 


98.5 


-0.7 


169.4(153.2 


± 18.3 


152.2 


147.1 


156.7 


155.3 


159.8 


-0.6 


208.7(192.5 


± 17.1 


201.1 


200.4 


205.1 


204.9 


218.3 


-0.5 


267.3(247.8 


± 19.3 


250.7 


247.7 


254.3 


247.7 


260.9 


-0.4 


320.3(296.7 


± 21.9 


296.0 


295.8 


294.4 


295.6 


293.0 


-0.3 


362.4(334.0 


± 22.0 


335.2 


338.6 


340.9 


329.8 


339.5 


-0.2 


399.1(366.7 


± 23.3 


368.5 


373.8 


369.3 


350.2 


369.1 


-0.1 


413.3(379.0 


± 25.9 


385.9 


389.7 


383.1 


369.5 


384.9 


0.0 


408.9(374.6 


± 23.0 


395.0 


407.4 


403.0 


388.0 


399.5 


0.1 


413.4(376.7 


± 23.1 


398.3 


409.6 


402.6 


372.2 


405.5 


0.2 


379.9(344.0 


± 25.8 


387.7 


398.3 


385.3 


354.9 


394.8 


0.3 


339.6(308.6 


± 21.1 


362.0 


368.7 


361.9 


329.8 


373.0 


0.4 


312.7(279.8 


± 21.5 


329.4 


331.6 


326.5 


296.9 


332.7 


0.5 


253.7(226.8 


± 18.4 


289.4 


286.8 


295.2 


254.7 


293.4 


0.6 


180.8(161.2 


± 16.3 


245.1 


239.0 


250.0 


206.3 


245.0 


0.7 


133.5(120.6 


± 12.6 


191.8 


193.4 


186.2 


159.8 


194.0 


0.8 


94.3(85.6 


± 12.8 


140.3 


139.9 


132.3 


110.2 


135.8 


0.9 


38.6(33.0 


± 13.4 


89.0 


90.1 


74.6 


57.4 


85.5 


1.0 


-19.7(-15.6 


± 14.8 


35.8 


37.4 


32.5 


3.0 


19.9 


1.1 


-32.6(-33.8 


± 14.4 


-11.9 


-19.9 


-26.5 


-39.5 


-29.8 


1.2 


-87.8(-80.7 


± 13.6 


-53.9 


-61.9 


-68.7 


-73.4 


-65.4 


1.3 


-117.6(-10S.8 


± 11.6 


-87.4 


-93.8 


-103.0 


-107.6 


-101.6 


1.4 


-135.0(-124.8 


± 13.3 


-113.9 


-121.6 


-127.5 


-129.6 


-121.1 


1.5 


-139.7(-125.3 


± 14.3 


-137.4 


-141.7 


-146.0 


-148.9 


-143.2 


1.6 


-134.0(-122.1 


± 14.3 


-152.9 


-153.7 


-163.7 


-160.9 


-159.4 


1.7 


-135.6(-121.3 


± 12.4 


-162.8 


-156.5 


-172.8 


-170.1 


-167.4 


1.8 


-140.0(-124.8 


± 13.6 


-165.2 


-158.8 


-168.3 


-169.9 


-170.7 


1.9 


-143.5(-127.3 


± 11.1 


-161.0 


-158.6 


-169.3 


-168.3 


-166.4 


2.0 


-135.7(-119.5 


± 12.5 


-156.8 


-155.7 


-160.2 


-158.2 


-158.0 


2.1 


-127.6(-111.7 


± 11.8 


-146.5 


-147.7 


-149.1 


-149.7 


-144.8 


2.2 


-116.8(-100.0 


± 11.3 


-134.4 


-133.7 


-141.9 


-132.7 


-137.0 


2.3 


-93.0(-Sl.l 


± 11.3 


-121.5 


-116.0 


-129.4 


-118.3 


-122.8 


2.4 


-84.4(-71. 


J) ± 8.6 


-110.9 


-100.8 


-114.4 


-104.1 


-108.8 


2.5 


-73.7(-63. 


3) it 7.2 


-97.4 


-91.1 


-96.9 


-91.7 


-91.3 


2.6 


-64.4{-54. 


7) ±8.6 


-83.6 


-77.1 


-81.0 


-77.3 


-79.8 


2.7 


-54.6{-47. 


2) ± 8.4 


-71.2 


-62.5 


-68.5 


-63.7 


-66.6 


2.8 


-38.1(-31. 


3) ± 8.1 


-61.3 


-51.0 


-55.9 


-52.4 


-54.3 


2.9 


-34.5(-29. 


3) ±7.1 


-49.2 


-43.4 


-44.0 


-44.2 


-44.8 


3.0 


-23.4{-19. 


J) ± 6.0 


-39.4 


-37.4 


-37.0 


-34.9 


-37.2 



Note. Genus values at a given threshold level of the SDSS galaxies in the BEST sample and five sets of the mock galaxies. 
The density fields were smoothed smoothing scale Rg is 6.0/i^^ Mpc. The genus values of the SDSS galaxies are systematic 
bias-corrected and the observed values are given in parentheses. 



